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^ ' Abstract. We consider a so-called random obstacle model for the mo- 

00 ' lion of a hypersurface through a field of random obstacles, driven by 

CO ■ a constant driving field. The resulting semi-linear parabolic PDF with 

' random coefficients does not admit a global nonnegative stationary so- 

lution, which implies that an interface that was flat originally cannot 
. get stationary. The absence of global stationary solutions is shown by 

' proving lower bounds on the growth of stationary solutions on large do- 

7—i , mains with Dirichlet boundary conditions. Difficulties arise because the 

• random lower order part of the equation cannot be bounded uniformly. 

> : 

^ ■ 1. Introduction 

H I 

I We are interested in the behavior of a moving interface F in a random 

medium, where F is a graph, i.e. defined as 

(1) r{t):={{x,y)eR^:y = u{x,t)} 

and the function u evolves according to the following equation: 

(2) — = u^^ix,t) + f{x,u{x,t)) + F in E x E+, 

(3) u{x,0)=0 
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where / G C^{M? x 0) is a random field representing the random medium 
and will be defined more precisely later on. Note that / is not restricted to 
be either positive or negative. F is a positive constant called "driving field." 
The objective is to prove that the solution of (2)-(3) does not get pinned, 
i.e. does not converge to a nonnegative stationary solution if F is above 
a critical value Fc- To this end, we will show that nonnegative stationary 
solutions on bounded intervals [—N, N] with Dirichlet boundary conditions 
get large with high probability as — )• oo. 

The main contribution of this paper is to show that a finite F is sufficient 
to keep the graph moving, even if it will have to pass through regions where 
f{x,u,uj) <C —1, provided the probability of finding such a region is small. 
As / can become arbitrarily big, one cannot find a deterministic subsolution 
that keeps moving, and instead probabilistic arguments are needed. 

The interest in the model stems from the theoretical analysis of the ef- 
fective behavior on large scales of models for interface evolution specified 
at a microscopic scale, which is at the heart of many problems in physics 
and material science. Of particular interest is the infiuence of material het- 
erogeneities, which are generally assumed to be random. Mathematically, 
this leads to studying the limit of evolution equations with rapidly varying 
random coefficients. In the case of dissipative equations, on which we fo- 
cus here, the randomness leads to new and interesting effects absent in the 
case of periodic coefficients, e.g. pinning and de-pinning for obstacles with 
a strength that cannot be bounded uniformly. If the strong obstacles are 
sufficiently rare, than the interaction through the Laplacian helps the graph 
overcome them although the total forcing f{x,u) + F remains negative near 
the obstacle. 

One example we have in mind as motivation are driven elastic systems, 
for a review of the research in physics and its possible applications we refer 
to [1]. For a survey of front evolutions in random media, with evolution laws 
different from the ones considered here, see e.g. the recent monograph [8]. 

The model (2) is obviously a gradient fiow for a random energy. In fact, 
it approximates a more geometric interface evolution law: 

In fact, if the hypersurface S is the boundary of the set Aj] then we can 
define for any bounded D QM? the energy 



where X € and denotes the 1-dimensional Hausdorff measure. 

Requiring that the first variation of that functional (with respect to inner 
variations, i.e. deforming the interface with the flow of a smooth vector field) 
is proportional to the normal velocity of the interface leads to forced mean 
curvature flow. 




V = K + f{X), 
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where k denotes the mean curvature of the interface (trace of the second 
fundamental form) and the scalar V is the velocity of the interface in the di- 
rection of the inner normal. This geometric evolution law leads to nonlinear 
degenerate parabolic equations, hence questions concerning the large-scale 
behaviour of solutions are related to homogenising such equations with pe- 
riodic or random coefficients. This is an active field of research (see e.g. [2], 
[6]) but many difficult problems remain open. Here we consider a modified 
evolution law: 

If we suppose that the interface is a graph is "fiat" (no overhangs, small 
gradients) then we can consider a semi-linear equation as in (2) as heuristic 
approximation of the evolution by forced mean curvature flow. 

This model, here called random obstacle model (ROM) because of the 
precise nature of the random nonlinearity f{x,u,io) used in this paper, is 
a special case of a class of equations sometimes called quenched Edwards- 
Wilkinson model which, for some choices of the random nonlinearity, is used 
in physics as a model for overdamped interface evolution in a random envi- 
ronment when "overhangs" can be neglected. For further comments on phys- 
ical properties and justifications of the model we refer to [1]. In particular, 
one expects that solutions move with a deterministic effective (large-scale) 
velocity for F larger than a critical forcing . For F slightly larger than F^ , 
the relation between the effective velocity and F — F^, is expected to be a 
power law. (See also [4] for the periodic case.). 

While there are important differences between the forced mean curvature 
flow and the semi-linear model (e.g. forced mean curvature flow can "wrap 
around" strong obstacles), we expect that the techniques we will develop 
when studying (2) will prove helpful in investigating more general models 
for interface evolution. This strategy was successful in the periodic case, 
where first the semi- linear case was solved ([4]) and then the results could 
be extended to graphs evolving by forced mean curvature fiow ([3]). 

One more reason why such models are of mathematical interest is the 
relation with "singular" homogenization problems, i.e. problems where the 
e-equation is of second order (possibly degenerate) and the homogenized 
equation of first order. Note that the effective velocity c(r/) of an interface 
evolving with average slope rj can be found by considering 

du 

— = u^xix, t) + f{x, r]-x + u{x, t)) + F, 
i.e. this can be seen as the "cell problem" for 

dv{y,T,u}) \ , ff -1 -If ^ ^ , Z7 
= (^Vyy{y,T,u}) + f{e y,e v{y,T,uj),uj) + F 

with T = e~^t, y = e~^x. 

The paper is organised as follows. In Section 2 we define the random 
obstacle model precisely and state our main results. 

In Section 3, we introduce an auxiliary model which is more suitable for 
explicit estimates and whose solutions can be related to solutions of the 
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original equation (2) by the comparison principle for parabolic equation. 
This auxiliary problem has the property that any of its stationary solutions 
u solve Uxx = —F away from the obstacles and is a convex function on the 
obstacles. This fact allows us to define a discretization, using that each 
solution is determined by its values when entering and leaving an obstacle. 
This yields a discretised path : Z — > (5Z characterizing each stationary 
solution. 

In section 4, we estimate the discrete Laplacian of v^{i) against the ob- 
stacles that sit above and below i ^ Tj and are approached by the path, 
i.e. lS.^v(i) + F < C£j [^5(-j)j(ti;) where F is a constant which can be chosen 
arbitrarily large. A technical problem is posed by the fact that the path may 
pass more than one obstacle above the same integer. 

In section 5 we estimate the probability of a discrete path being "com- 
patible" with the random environment. This probability can be estimated 
against an auxiliary random measure on paths: 

P ({(J : u{uj) compatible with < C^^P [{Aav^ {i)}^J-_\^^^ , 

where Z is a normalization (corresponding to the partition function in sta- 
tistical mechanics). 

In section 6 we conclude that the probability of a nonnegative solution of 
the Dirichlet problem to cross KN — K\x\ is 0{e~'-^^). The key observation 
is that for such a path N~'^ Yli^~N+i i^dV^ii) + F^ must be large, which is 
very unlikely under the auxiliary (product) probability measure. 

Finally, we show by invoking the comparison principle for semi-linear par- 
abolic equations that these results for large N imply non-existence of global 
nonnegative stationary solutions. This implies that for a solution u of (2), 
(3) and all x G M it holds that liuit^oo u{t, x,uj) = +oo almost surely in uj, 
i.e. the interface cannot be stopped by the obstacles. 

Acknowledgements The second named author would like to thank Enza 
Orlandi and Michael Scheutzow for helpful discussions. The authors ac- 
knowledge gratefully the hospitality of the Max Planck Institute for Math- 
ematics in the Sciences Leipzig. 

2. Results and Definitions 

2.1. The random field /. Here, the field / is negative on "obstacles" in 
which are random in strength, but positioned on a lattice. More precisely, 
we make the following assumption: 

Definition 2.1 (Obstacles). 

(1) Let Z* := Z + 1/2. We assume that the obstacles lie on a lattice 
£ := Z X Z* where for convenience (^ijO^gg ^gg. denotes the nodes of 
this lattice, i.e bij := (i, j). 
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Figure 1. The obstacles 



(2) Let 5 < 1/2 and define Q5 (0,0) := [-6,6]'^ , and Qs{i, j) := Qs{0,0)+ 
bij. Then the obstacles, i.e. regions where / < is possible, are given 
by the Qsii^j), see also Figure 1 

In order to obtain existence and regularity of the solutions, the nonlinear- 
ity f{x,y) should be sufficiently regular, hence in order to define / we have 
to smooth out the obstacles. 

Definition 2.2 (Random field). Let (j) € be a nonnegative function such 
that its support is contained in cube Qs{0,0). 

Let (H«, j)(w))(i,j)gzxz* be a family of independent identically distributed 
exponential random variables, i.e. there exists Aq > such that for r > 

¥{l{i,j){uj) >r} = e-^°\ 

Let S be the set of the obstacles, i.e. S := U(ij)GZxZ* (Qsi^ij)^ , then 
the field / is defined the following way: 



f{x,s) = g{x,s) - ^ l{ij)(t){{x,s) -bij) 
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where g is a non-negative function chosen so that the field has mean zero in 
a suitable sense: 



g>0 in 
lim {2Lf / f{x,s)dxds = 

J\-L,LP 



Remark 1. (1) As = j, the law of large numbers implies that 

a possible choice of g is 

g{x,s)= j(f>{{x,s) -bij). 

(i,j)GZxZ* 

(2) The results on non-existence of nonnegative stationary solutions hold 
for any i.i.d. random variables such that there exists Aq > 
with 

F{l{i,j){uj) >r} < e"^'"'. 

(3) As we are only interested in the combined effect of f{x,s) and the 
constant forcing F, the mean zero property of the random nonlinear- 
ity is just a normalisation. 

(4) In our analysis, the shape of the obstacles {supp{(j))) plays no role 
and the results will stand as well if we consider a random field like 
e.g. 

f = g{x,s)- ^ l{i,j)(l)ij{{x,s)) 
(ij)ezxz* 

where 4>ij are smooth functions uniformly bounded and such that 

supp{4>ij) C Qs{bi,j) ■ 
2.2. Results: We consider the stationary version of (2) with Dirichlet 
boundary conditions: 

(4) u^^ + f{x,u,u}) + F = in [-N + 5,N - 5] 

(5) u{-N + 6) =u{N -6) = 

Theorem 2.3. Let u{uj) solve (4, 5). Then there exist Fq > 0, C and K 
such that for F > Fq and N sufficiently large 

F{{uj\u{x,io) > {K{N -1) - K\x\)+ on [-N + 6,N -6]}) > l-Ce'^, 
where denotes the positive part of a real number a. 

Corollary 1. Let F > Fq, with Fq as in Theorem 2.3. 

(1) There is almost surely no global nonnegative stationary solution of 
(2). 

(2) Let u solve (2), (3). Then 

lim u{t,x,uj) = +00 for all x S M 



holds with probability one. 
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3. Blocked path and auxiliary problem 

In this section we define a auxiliary problem that we will constantly use 
along this paper. We will denote by xb the characteristic function of the set 
B. 

Definition 3.1 (Auxiliary field). Let 

A:=W^\{\J{i-5,i + 5) xR} 

■.= ^^ \ - 5 - e,i + 5 + e) x R} 

and define 

7(x,s):=- ^ l{i,j)(l){{x,s) -bij). 
{i,j)ez*xz* 

Let us now consider the following auxiliary problem 

dv ~ 

(6) = + f{x, vit, x)) + Fx'a{x) 

(7) v{0,x)=0, 

where Xa ^ smooth function such that XA^ 1^ Xa — Xa- e is a small 
parameter which will be fix later on. 

To visualize the new random field defined by g = / + Fxa{x) see figure 
2. Note that it is differentiable in x and s. _ 

Observe that, as the obstacles are negative, / + Fxa ^ f + F- Therefore 
the comparison principle for the parabolic equation (see section 3.1) implies 
that solutions of the auxiliary problem remain below solutions of the orig- 
inal problem. Hence existence of a nonnegative stationary solution for the 
original problem implies existence of one for the auxiliary problem. By con- 
traposition, nonexistence for the auxiliary problem implies nonexistence for 
the original problem. 

Stationary sub/supersolutions can be constructed as piecewise quadratic 
functions. For any F we can construct the graph of such a solution (also 
called "paths" to emphasize the analogy with a stochastic process). 

Definition 3.2 (blocked path). A graph {x,v{x)) is called blocked path if 
and only if t; G Cl^jM), and 

(8) Vxx = -Fxa{x) in {i + 5,i + l-6), 

(9) Vxx=^^l{'i,i){'^)(t}i,j{x,v{x)) in {i-5,i + 5). 

where ^i^j{x,s) := (/>((x,s) - 6jj). 

Observe that the path for a; G (z + (5, i + 1 — 5) is uniquely determined by 
the boundary values v{i^S) and u(i + l — (5), because it solves a /inear elliptic 
equation there. But note that, for a given realisation of the random field. 
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Figure 2. Mapping of the obstacles for the auxihary problem 

there may be more than one blocked path, as equations like u^x = fix,u) 
do not have unique solutions without further conditions on the nonlinearity. 

Remark 2. /,From Definition 3.2, we see that f is a convex function in 
{i — 6,i + S) and hence we have 

v{i + 6)> v{i -6) + 26v'{i - 6) 

Let us now define some discrete quantities that we will use throughout 
the paper. 

Definition 3.3. Let v{i) and v^[i] be defined as follows: 

v{i) := v{i — 6) + 25vx{i — 5)-, 



inf{j G > v{i) - -} & S'i 



We will need the following Lemma. 

Lemma 3.4. Let v be as in Definition 3.2 and v, he in Definition 3.3. 
Denote by the piecewise linear interpolation ofv^, and by w the piecewise 
linear interpolation of v. Then v + 5/2 > iju^ , and v > w. 
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Proof. First, note that convexity of v in [i—6, i+6] implies that v{i) < v{i+5). 
Let 

Ii:= {i-l + 6,i + 5) 
and let the auxiliary function w be the solution of 

Aw = -Fl[i_i_^_s^i_s] on li 

w{i — 1 + 5) = v{i — 1 + 5), w{i + 5) = v{i). 

This function is on its domain and solves the ODE w^x = —F on (i — 1 + 
S,i — 5). (Here x is considered as "time"). Suppose w^i — 6) > v{i — 5). Then 
Wx{i — 6) < Vx{i — 6), and integrating the ODE backwards in x we obtain 
w{i — 1 + 5) > v{i — 1 + 6), a contradiction. Assuming w^i — 6) < v{i — 5) we 
obtain a contradiction in a similar way, and we conclude w{i — 5) = v{i — 6). 
This implies w = v on [i — 1 + 6,i — 5] and (by convexity ol v on [i — 6,i + 6]) 
w < V on [i — 1 + 5,i — 6]. 
Now consider 

Aw = on /j 

w{i — 1 + 6) = v{i — 1), w{i + 5) = v{i) 

Clearly w is the piecewise linear interpolation of v. 

As Aw < Aw and w < w on dli, the comparison principle for the Laplace 
operator implies w > w, so v > w > w. The conclusion for w^ follows 
immediately. 

□ 

3.1. Existence and uniqueness for parabolic equations. 

Lemma 3.5. There exists a global classical solution of the parabolic Cauchy 
problems (2), and (6) with initial conditions which are uniformly bounded 
and locally C^. The solutions are unique. If < vq < uq, v solves (6) with 
initial condition uq, u solves (2) with initial condition uq, then v < u. 

Proof: For M G N, replace l{i,j){uj) by l^^{i,j) := M A l{i,j), where A 
denotes the operation a A b := inf{a, 6}.The corresponding fields f^^, 
are uniformly bounded and uniformly Lipschitz in s. Therefore we can apply 
the Banach fixed point theorem in L°° in order to obtain a local in time 
solution, which, by local parabolic regularity, is classical. It can be extended 
as the nonlinearity is uniformly bounded. Hence a global solution u^{x,t) 
exists. Note that by the comparison principle u^^ is a positive monotonic 
non-increasing function of M i.e. n*^ > > for > M , so u{x,t) : = 
limM^ooU^ ix,t) exists. Applying regularity locally, (where the obstacles 
are bounded) we obtain that the limit is a classical solution. 

4. a priori ESTIMATES ON v{i) AND v^[i] 

In this section, we establish some a-priori estimates on v{i) and v^[i]. 

First we show a lemma which allows to estimate the discrete Laplacian of 
"0 at i (which involves i,i + 1 and i — 1) by something that depends only on 
the obstacles above i. 
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Lemma 4.1. Let v{i) defined as in the previous section, and define the 
discrete Laplacian as 

Adv{i) := v{i + 1) - 2v{i) + v{i - 1) = {v{i + 1) - v{i)) - {v{i) - v{i - 1)) 
Then 

-26[v,,ii-l + 6)-v^ii-l-6)] < Adv{i) + F < {l + 26)[v^ii + 6)-v^{i-6)]. 

where F{1 - 2{6 + e)) < F < {I - 26)F for the e > in Def. 3.1. 

Note that our discretization, using the tangents, imphes that the discrete 
Laplacian does not necessarily satisfy the same lower bound as the Laplacian 
of the original path. 
Proof: 

Step One : Upper Bound 

As a preparation, let us recall some formulas satisfied by v. 
Since v satisfies (8), we have for alii G Z 

(10) 

Va:{i + 1-S) -Va:{i + S) = -F XA{x)dx 

J i+5 

(11) 

v{i + l-5)-v{i + 6) = {l-25)v^{i + 5)- F f ([ XAix)dx] ds 

Ji+5 \Ji+5 J 

Let us define 

i-i+l-S 

F:=F XAi^)dx. 

J i+S 

Observe that since Xa(^ +p) = Xa(^) integer p, F is independant of 

i G Z. Moreover 

F{1 - 2(5 + e))<F <{1- 26)F 

since xaA^) < Xa(^) ^ Xa{x). 

Using now (10), the definition of v{i + 1) and (11) we see that 

ri+l—5 / l-s \ 

v[iJ^l) = v[i + 5)+v^{i + 6)-F / XA{x)dx] ds + 25{v^{i + ls) -v^{i + 5)) 

Ji+5 \Ji+5 J 

ri+1—5 / l-s \ 

= v{i + 6)+ vS + 6)-F i XAi^)dx ds - 26F. 

Ji+5 \Ji+5 J 

Therefore, 
(12) 

v{i+l)-v{i) = v{i+6)+Vx{i+d)-F i XA{^)dx] ds-25F-v{i). 



I i+5 \Ji+5 

Observe that since Xa(^ +p) = Xa(^) fo^^ integer p we have 

F i XA(.x)dx] ds+26F = F / I / XA(.x)dx) ds+26F. 

Ji+5 \Ji+5 ) Ji-1+5 \Ji-l+5 J 
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Hence, from the definition of the discrete laplacian and using (12) it follows 
that 

(13) Adv{i) = v{i + S) + v.^{i + 6)-v{i)-v{i-l + 6)-v.j:{i-l + 6) + v{i-l) 

Using now the definition of v{i) and the convexity of v in (i — 5, i + 5) for 
all i € Z we see that 

v{i + 5)+ v^{i + 5)- v{i) < v^{i + 6) + 25iv^{i + 5)- ^^^(i - 5)) 

-v{i-l + S) + v{i-l) <0. 

Hence, 

Adv{i) < (1 + + S)- 26v^{i -6)- v^{i -1 + 6). 

Using now (10) it follows that 

Advii) < (1 + 2<5)(z;^(i + 5) - v,{i - 5)) - F. 

Step two: Lower bound 

From formula (13) we have 

(14) Adv{i) = v{i + 6)-v{i) + Vx{i + 5)-v{i-l + 6) + v{i-l)-Vr,{i-l + 5) 

Since v is convex in {i — 5,i + 5), we have v{i + 5)—v{i) > and Vx{i + S) > 
Vx{i — <5)- Therefore we have 

(15) Adv{i) > Vx{i -6)- Vx{i - 1 + 6) - v{i - I + 6) + v{i - I). 

Using now (10), the convexity of v in (i — 1 — 5, i — 1 + 6) and the definition 
of v{i — 1) it follows that 

(16) Adv{i) >-F- 25[vx{i -l + 6)-Vx{i-l-S)]. 

Now we proceed to estimate the change of the discrete gradient 

k{i) := Vx{i + 5) - Vx{i - 5) 

in terms of the obstacle strengths above i. Observe that always A; > by 
convexity. If the gradients are very steep, the path will pass through several 
obstacles above the interval [i — 6,i + 5\. The number of obstacles passed and 
the time spent in each of them (i.e. the Lebesgue measure of its image under 
the inverse mapping) can be estimated in terms of v'{i — 5) and v'{i + 5). 

Lemma 4.2. Let v he a blocked path, i ^"L and and assume that k{i) > 0. 
Set M := sup{\vx{i — \vx{i + S)\} then we have 

v[i)-ASM<j<v{i)+A&M 

Proof: 

Step 1: As f is convex on [i — + 5], the gradient is monotone, hence 
|^^x(f^)| < M for all X G := [i — (5, i + 5\. As a consequence, we have on 

m 

v{i) - 6M < v{x) < v{i) + 6M, 
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As \v{i) -v{i- 5)\ < 26M, \v{i) - v{i - 6)\ < 5M, we obtain 

\v{x) - v{i)\ < 46M on[i-6,i + 6]. 

Step 2. Define the time spent by the path in the j-th. obstacle above i as 

S-j ■.= \{x: v{x) e [j-5,j + d]}\, 

where \ A\ denotes the Lebesgue measure of the set A and j € = 1/2 + Z. 
Note that by convexity Vx changes sign at most once, hence each Sj is the 
union of at most two intervals, moreover Sj = ^ if \j — v{i)\ > 46 AI Hence, 
as for X € Vxxi^) < lihj) on obstacle j and zero else, 

Vx{i + 5) - Vx{i - 5) < ^ l{iJ)Sj. 

v{i)-i5M<j<v{i)+i5M 

where j G Z* 

Step 3. Note that k < 2M. As the gradient is monotone on /(i), there 
exists a f such that \vx{f)\ = M - k/3 and \vx{x)\ > M - k/3> M/3 > 
on where 



m 



[f,i + 6] if M = \vxii + S)\ 
[i-6,f] if M = \vx{i - 6)\. 



As the gradient does not change sign on the sets Sj := Sj n are 
intervals. Moreover, 

I c I < = ^ 

' ^' - M/3 M 

as \vx\ > M/3 on Hence 

v{i)-4SM<j<v{i)+4:5M v {i)-4:5M<j<v(i)+45M 

and the result follows. 

□ 



Remark 3. Note that in the case where k(i) > 1 then the corresponding 
M{i) > ^. Indeed, by definition of M{i) and k{i) we have 2M{i) > \vx{i — 
5)1 + \vx{i + -5)1 > Vxli + 6)- Vxii -5) = k{i) > 1, i.e. M{i) > 1/2. 

Combining now Lemmas 4.1 and 4.2 we deduce the following estimates, 
which allow to estimate the discrete Laplacian of the blocking path (7;^[j])jg[_jv,Ar]nz 
at a site i against a normalized sum of random variables. 
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Lemma 4.3. Let v be a blocked path. Then for all i € [— + 5, N — 5] (IZ 

there exists M{i),M{i — 1) > ^ such that following holds 



Advii)+F< {1 + 26) 



> -26 



360(52 



25(4M(i) + i; 



360(5^ 



E 



v{i)-&{AM{i)+\)<j<vii)+5{AM{i)+\) 



E 



2<5(4M(i - 1) + i) 

^ ^ ' ^' v(i-i)-S{AM{i-l)+\)<j<v{i-l)+5{4:M{i-l)+\) 



l{i-l,j){i0) 



where F := F - {I + 26). 
Proof: 

Let us first start with the proof of the upper bound. Observe first that 
v{i)-^-<v'[i]<v{^) + ^-, 

which implies that 

l^dv^i] -26< /\dv\i\ < /:^dv^[i\ + 26. 
Therefore using Lemma 4.1 we have 
(17) AdV^[i] < {l + 26)kii) - F + 26. 

with k{i) > 0. By Lemma 4.2 and Remark 3, for k{i) > 1 there exists 
M{i) > ^ so that 



k{i) < 



186 
M{i) 



E 



l{i,j){u). 



v{i)-45M{i)<j<v{i)+i5M{i) 

So we easily see that 

(18) < ^^^''^^"<'> + f > V 

^ ^ ^ ^ - M(i)(4M(i) + i)<5 , 1 ^ . 



l{i,j){uj). 



Therefore, since M(i) > ^ we have 
180(52 



(19) A;(i) < 



J] /(i,i)(a;). 

2^" t)'5[j]-(4Af(i)+i)5<j<i}«[i]+(4Af(i)+i)(5 



(4M(f) + i)(5 



Hence, for all k{i) > 0, we have 
180^2 

k{i) < 1 + 



{AM{i) + i)(5 



l{i,j){uj). 



2^ i)'5[j]_(4M(j)+i)5<j<S'5[i]+(4Af(i)+§)<5 

and the estimate follows . The lower bound is treated in a similar way. 



□ 
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5. Probabilistic Estimates 

We first recall a standard fact for the Laplace transform of independent 
exponential random variables and random variables with distribution func- 
tion bounded by an exponential. 

Lemma 5.1. (1) Let {XjjjgTv independent identically distributed ran- 
dom variables such that for a parameter Aq and a constant C > 

lP[^o >r]< Ce-^«^ 
Then we have for any A < Aq and L £ N, L > 2, 

Ao 



(20) 



(21) 
(22) 
(23) E 



E 



E e^^i 



< C 



An -A 



< c 



An 



Ao-A 

41n(4/3)A 
3Ao 



for L > 2, A G (2/3Ao, Ao) 



(2) Let {Xjjjgjv be independent exponential random variables with pa- 
rameter Ao > 0. Then (21)-(22) hold as equalities with C = 1, while 
(23) holds as inequality with C = 1. 

Proof: 

We first show 2. The first equality is standard, the second follows by using 
independence. For the third, note that by concavity of ln(l — x) on [0,3/4] 



ln(l — x) > — xln(3/4) for x G 
3 



Using independence and this concavity estimate with x = Xq/{XL) 



E 



A, 







An 



A 

L . 



4A 
3Ao 



In order to show 2., it is sufficient to prove the first inequality, the others 
then follow as in the previous case. For (21) note that the expectation of 
a random variable is the Riemann-Stieltjes integral with the distribution 
function as integrator. Now integrate by parts and use that the integrand 



gAx monotone. 



□ 



Remark 4. Observe that the above estimate on the Laplace transform of 
Sl is independent of L. 



Let us define Sm by 



SM{uj){i,j) :-- 



The we have the following Corollary: 



'M<j-l<M 



RANDOM COEFFICIENTS 
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Corollary 2. For any discrete function j{i) : Z ^ Z, the random variables 
{S'M(w)(i, j(^))}«ez ore independent and identically distributed. Moreover, 
there exist constants C, A which depend only on Aq such that 

F (SM{i^){iJ{i)) > <e^-^^ 

Proof: 

The first assertion is obvious. The second is a consequence of (21) and (23) 
and the exponential Chebyshev inequality with a parameter A S (2/3Ao, Aq). 

□ 

Let us now estimate the probability of a blocked path with boundary 
conditions on [— A^, A^] to be compatible with the l{i,j). 

Definition 5.2. (blocked Dirichlet path) Let v{-N + 5) = v{N - 5) = 0. 

Moreover, let v solve (8) for — < i < N — 1, and let v solve (9) for 

-N + l<i<N-l. 

Extend v to [-N -6,N + 6] by 

v(^x)=v'{-N + 6){x + N-6) on [-N -6,-N + 6] 

and 

v{x) = v'{N -6){x - N + 6) on [N-6,N + 6]. 

Remark 5. 

(1) Note that this path solves (9) for -N < i < N ii we set l{i,j) = 
for i = —N or i = N. 

(2) If u > on [-N + 5,N-6], then 

0>v{x)>-26FN for x e [-N - S, -N + S] U [N - 5, N + 5]. 

Definition 5.3. Let : [-N, N]nZ ^ SZhe a discrete path. We call the 
path compatible with a random obstacle configuration if there exists a (not 
necessarily unique) path as in Definition 5.2 which is mapped to under 
the discretization defined in Def. 3.3. 

Note that the discrete path is fixed. Whether it is compatible or not 
depends on the configuration of the random field. 

Lemma 5.4. Let (17, J-", "P) be a probability space and let l{i,j){uj) be i.i.d. 
exponential random variables with parameter Aq > and let be a discrete 
path with fixed boundary conditions 

v\-N + 6)=0, v\N) = b for some b G [-FN, FN]. 

Then there exist constants C{5,Xq), Ai((5, Aq) independent ofb such that we 
have for F sufficiently large 

F[v^ compatible, v^{N) = b] := ¥b[v^ compatible] < e^^e"^i^^~+i 
with F as in Lemma 4-3 
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The previous estimates bounds the probability of the random obstacle 
configurations such that a fixed discrete path is compatible with the random 
environment. In order to prove that the probability that there exists some 
compatible nonnegative path is small, we would have to sum over all possible 
paths, each weighted with the right hand side of the previous estimate. It is 
complicated to bound these sums, because the number of possible discrete 
paths grows faster than exponentially in A^. Fortunately, most of them are 
extremely unlikely to be compatible. In order to quantify this, we define an 
auxiliary probability measure on discrete paths. 

Definition 5.5. 



oo 




g-Ai|<5fc+F| 



A;=— oo 



The normalisation constant is obtained by summing over all possible dis- 
crete paths for fixed boundary conditions. This is equivalent to summing 
over all discrete Laplacians. Note that Z is bounded from above and below 
by constants independent of F. 

Note that the law of the positive and the negative part of AdV^{i) + F 
under P is that of (discretized) independent exponential random variables. 
In particular, probabilities of sums of the discrete Laplacians have certain 
exponential moments and can be estimated by large deviation techniques. 

Corollary 3. With P as in Def. 5.5, there exists A'^o(^O)^) such that 

¥b[v^ compatible] < e^^FlAdV^] 

for N > No. 

Proof of Corollary 3: 

We suppose that Lemma 5.4 holds. Then 

ni^'' compatible] < e^^e"^^ l^'^^'W+^l 

= e^^ {ZY Z-'^^e~^^ ^-^+1 |A.f)^«+F| < ^NC^^ 
Z 

for sufficiently large. Here we can choose e.g. 

C = 2C + 2ln{Z). 

Proof of Lemma 5.4: 

In order to simplify notation we write 

S^s{uj){i) := SM{v^){^){i,v^{'i'))- 
We write the absolute value as sum of positive and negative part. 
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By Lemma 4.3 we get that there exist universal positive constants Co such 
that the fixed discrete path is compatible only if 



N-2 



^ G n {A^s^+{i) n A,s^_{i)) n A^s^+{N - 1) n A^s^_{-N 

A,s^+ii) := |w: Co(^Adv\i) + F)_^<S^siu;){i) 
A,s^_ii) := L: Co[Adv'ii + l) + Fy<S,s{oj)ii) 



B„s{i) := {A^s^+{i)nA^s^_{i)). 
Note that 

B,s{i) C ls,s{i) > % (A,v'{i + 1) + f) + % ( Arft;''(i) + F 



and we estimate with the help of Corollary 2 for i G {— A^ + 1,...,A^ — 2} 
P(S,.(i)) < /-If ((A.^)*W+i^)^ + (A..^(i+l)+F)_) 



for constants C and Ai depending only on Aq but not on F. 
Moreover, for i = — 1 we obtain 

and for i = —N + 1 we obtain 

P(A,.,„(-iV+ 1)) < e^-^(^^^^(-^+^)+^)- 
The events B^s{i) are independent for different i, hence 

Pb[?;'^ compatible] < F{A^,s ^^{-N + l))F{A^,s ^^{N - 1)) ]J F{B^s{i)) 

i=-N+l 

The claim follows now by choosing Ai = ^r^. 

Remark 6. Note that the 1-1-correspondence between second derivatives 
and paths with Dirichlet boundary conditions allows us to express each path 
uniquely through its discrete Laplacians and thus estimate its probability 
with the help of the previous lemma. 

As a consequence the discrete Laplacians on average much larger than —F 
are extremely unlikely. We will show that nonnegative paths that cross the 
"triangle KN — K\x\ require such unlikely values of the discrete Laplacian. 



□ 
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6. Final Argumentation 

6.1. Some formulas on discrete path and Comparison of two paths. 

In this section, we recall some well known formulas for discrete paths and 
their discrete derivatives. The proofs are straightforward computations and 
therefore omitted. 

Let us first recall some basic formulas satisfied by a discrete path z defined 
in Z X M. 

Lemma 6.1. Let us denote V^z[£ + 1] := z[i + 1] - z[i] and V z\i + 1] := 
z{i + 1] - z{i + 2]. Then fori€Z we have 



Vz[l + 1] = ArfzM + Vz[l\ = Adz[i] + V'z[l] 

i=l 

V^z[i + 1] = Adzli] + V^zli] = ^ Adz[^ + Vh[k]. 



i=k 



Ai + 1] - ^[0] = Y.Y. ^rf^C-?'] + + 
i=\ j=l 

e+1 I i 

,{l + 1] - z\k\ = (zH - z\i -\\)=Y E ^'^^[J'] + + 1 - + !]• 

j=fc+l j=fc+l j=fc+l 

(Hi) 

t 

V^z[0] = A,z[l] + V'z[l] = Y,^dz\i] + Vz[t], 

i=l 

I 

V^z[A;] = Arfz[A: + 1] + Vz[A; + 1] = ^ Adz[i] + V^z[£]. 

i=k+l 

(iv) 

e-i I 

z[o] -z[i+i] = Y,Yl + + 

i=0 j=i+l 

I e-i e 

z[k] - z[l + 1] = - z[i + l]) = Y.Y. ^d^^^^ + + 1 - k)V'z[l]. 

i=k i=k j=i+l 

(v) 

V'z[^+ 1] = -Vz[i] 



Let us now define what we mean by "crossing. 
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Definition 6.2. Let zi and Z2 be two given paths in Z x R. We say that 
zi cross 22 if and only if there exists i € Z such that zi[i] > Z2[i] and 
+ 1] < Z2[i + l]. 

We will apply this to the discrete path and the triangle zxii) '■= NK — 
K\i\. 

6.2. Proof of Theorem 2.3. First we state a trivial fact for discrete sums. 

Lemma 6.3. Let aj be nonnegative numbers, then 

N N N N 

(24) EE"^ = E-?«^-^^E"^- 

i=l j=i i=l i=l 

We will show that paths that remain nonnegative but cross the triangle 
Zk require values of the average discrete Laplacian which are very unlikely 
under P. In order to do so, we distinguish cases: Either the path is above the 
triangle near one of the two endpoints of the interval [—N, N] and crosses at 
the interior, or it crosses at N or —N. In both cases, this implies information 
on the gradient. Note that the nonnegativity of the original subsolution 
does not imply the nonnegativity of the discretized path, but only that the 
discretized path is larger than —6FN, 6 times the minimal possible gradient. 
In particular, it implies that the terminal value b of the discretized path is 
in [-6FN,0]. 

Notation: As only discrete paths appear in the following estimates, we 
will write v[i] for v^[i] In order to simplify notation. 
If —V^'v[—N] < K, then by Lemma 6.1 

v[0]-v[-N]= E ^dvijj-NV'-vi-N]. 

i=-N+lj=-N+l 

Since f [— A^] = and rewriting the double sum the right way, it follows that 

-1 

-FN < v[0] <NK+ {-i){Adv\i]). 

i=-N+l 

After adding and subtracting F in each term in the summation 
-FN<NK+ Y {-i){^dvm + F) - F ^^^ ~ 

i=~N+l 

SO, invoking (24) it follows that 

F^^^^^^-{F + K)N<2{N-1) J2 {^dv[^+F)+. 

i=-N+l 

By definition of F, we have 

F> F{1- 2(5 + e))- (1 + 25). 
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Therefore for e small, says e < 6 and F such that F > 2 j^j^g^ we achieve 

- 2 

Whence 

F ^^^~ -{F + K)N < 2{N - 1) {^dv\i]+F)+. 

i=-N+l 

This implies that for N large and K fixed 

1 1-25 



2(N ^, 

As the (A(if[i] are independent random variables under the auxiliary 

probability measure P defined in Def. 5.5 which have exponential moments 
bounded as in (21), we can derive an upper bound for the large deviations 
principle: (For the basic form of the large deviations principle needed, see 
e.g [5] Ch. 5.11) Let 

where /i := Aq^ with Aq as in Lemma 5.1. (I.e. for exponential random 
variables /i is the expectation of (Arfu[i] + F)^ under P. Note that /i is 
decreasing in Aq.) Then, by the large deviations principle, for any r/ > 
there exists iVo ^ N such that for all N > Nq 

r J' (A.^[^l - P). >- ^ .-(-^(^)-). 

where C is the constant in the bound (21). (C = 1 for exponential random 
variables.) Now choose F sufficiently large such that 



CC-I 



e V B 7 < 1, 

where the constants are defined in Lemma 5.4. 

Then there exists a constant C3 depending on Aq and 5 such that for 
sufficiently large 

P(case 1) < 6"^=*^. 
The case V^v[N] > —K is done in a similar way. 

Second case: —V^'v[—N] > K,V^v[N] < —K. This implies that the path 
has to cross the triangle inside the interval [—N, N]. Suppose the path crosses 
zk on [— A^, 0], the other case is follows by symmetry. Then there exists 
A^i, -Af < A^i < 0, such that -V'vlNi] < K and v[Ni] < KN. Then by 
Lemma 6.1 

N-l i 

v[N]-v[Ni]= Y E ^dv[j] - {N - Ni)V''v[Ni], 

i=Ni+l j=Ni+l 
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SO 

-FN < v[N] < 2KN+KN+ ^ {AMj]+F)-F ^^ ~ ^'^^^^ ~ ~ ^\ 

i=Ni+l j=Ni+l 

which imphes 

_N(N-^^ ' - N-i _ 

F ^ ^ ) -iF+3K)N< ^ i^Mj]+F)<2{N-l) ^ (A,(^[i])+F)+, 

i=Ni+lj=Ni+l i=-N+l 

i.e. for N sufficiently large 

N-l 



(1 - 25)F 



2(N ^, 



Now we can repeat the probabilistic argument from the first case. 

Finally, we sum over all possible —FiV values of the the terminal condition 
h. This sum grows linearly in A^, hence using the exponential decay of the 
probabilities we obtain that there exists C/^ib^ Aq) and -Fo('^) ^^o) such that for 
F>Fo 

P(c;j : compatible and crosses zk) < e"*^*^. 
Now we conclude with Lemma 3.4. 

□ 

6.3. Proof of Corollary. Define as the solution of the initial-boundary 
value problem 

= v^,{x,t) + f{x,v''{x,t)) + F in (-N + 6,N-6), 

v^{-N,t) = u = v^{N,t) = 
v'^{x,0) = 0, 

and let u{x, t) solve 2. The comparison principle for parabolic equations 
implies that v'^{x,t) < u{x,t) for x G [-N -6,N + 6], t > 0. 

Moreover, v'^{x,t) v^nti^) as t ^ oo, where fi^at(^) ^ stationary 
solution of the Dirichlet problem. 

Note that dtv^{x,t) > as dtv'^{x,0) > 0, and the time derivative w := 
dtv'^ solves 

dtw = Aw + V{x)w, 

where the potential V{x) = ^{x,v^ {x,t)) is bounded on compact subsets 
of M^. (Note that a; is a fixed parameter here, < FN'^l^_]^j^f^,so only 
obstacles within [—N, N] x FN'^ can occur, but these are bounded for uj 
fixed.) 

Now a linear parabolic PDE with sufficiently regular potential V{x) and 
nonnegative initial condition remains nonnegative: w = e~*ll^ll°°-u; solves 

dtw = Aw + V{x)w, V <0 
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with initial condition 5; > 0. So the classical parabolic comparison principle 
([7]) implies w >0. 

By Thm. 2.3 and the first Borel-Cantelli Lemma (see e.g. [5]), 

P (a; : v^^^^{0) < KN for infinitely many N) = 0, 

so there exist almost surely arbitrarily large N such that 

hminf n(0,t) > lim t;^(0,t) = vg^,{Q) > KN, 

which implies 

(25) liminf u(0, i,a;) = +oo 

i— >oo 

with probability 1. By the comparison principle, this contradicts the exis- 
tence of a global nonnegative stationary solution. 

Moreover, by arguments as in Lemma 3.4, ((25)) holds for x G [—1, 1]. As 
the distribution of the obstacles is invariant under translations in x-direction, 
((25)) holds for x e M. 
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